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Abstract: We investigate the power-suppressed corrections to the fragmentation functions of 
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by selecting final-state particles in the current hemisphere in the Breit frame of reference. Our 
method is based on an analysis of one-loop Feynman graphs containing a massive gluon, which 
is equivalent to the evaluation of leading infrared renormalon contributions. We find that the 
leading corrections are proportional to 1/Q 2 , as in e + e~ annihilation, but their functional forms 
are different. We give quantitative estimates based on the hypothesis of universal low-energy 
behaviour of the strong coupling. 
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1 Introduction 



The study of final-state properties in deep inelastic lepton scattering (DIS) has received a great 
impetus from the increasing quantity and kinematic range of the HERA data. Amongst the 
most interesting quantities being studied are the fragmentation functions [1], which specify the 
single-hadron momentum distributions resulting from the fragmentation of the struck parton. 
Although these functions cannot be calculated using perturbative QCD, their asymptotic scaling 
violations (logarithmic Q 2 dependence) can be predicted and used to measure the strong coupling 
a s , in a similar way to the scaling violations in the deep inelastic structure functions. In addition, 
the fragmentation functions and their scaling violations can be compared with those measured 
in other processes, such as e + e~ annihilation. 

One problem with the measurement of a s using scaling violation, in either structure or frag- 
mentation functions, is that there is Q 2 dependence associated with power-suppressed (higher- 
twist) contributions, in addition to the dominant logarithmic dependence. These contributions 
need to be estimated in order to make use of the wide Q 2 coverage of HERA. 

Recently so-called 'renormalon' or 'dispersive' methods of estimating power-suppressed terms 
have been suggested. By looking at the behaviour of the QCD perturbation series in high orders, 
one can identify unsummable, factorially divergent sets of contributions (infrared renormalons 
[2, 3, 4]) which indicate that non-perturbative power-suppressed corrections must be included. 
The Q 2 -dependence of the leading correction to a given quantity can be inferred, and by making 
further universality assumptions one may also estimate its magnitude. Tests of these ideas 
provide information on the transition from the perturbative to the non-perturbative regime in 
QCD. In particular, one can investigate the possibility that an approximately universal low- 
energy form of the strong coupling may be a useful phenomenological concept [5]-[8]. 

Such an approach has been applied to a wide variety of observables, including DIS structure 
functions [8]-[12], e + e _ fragmentation functions [13, 14], and event shape variables in e + e~ 
annihilation [15]— [21] and DIS [22, 23]. Comparisons with experimental data [24]-[27] have 
been encouraging. 

Here we extend the same method to DIS fragmentation functions, considering in the present 
paper the contribution which is non-singlet with respect to the incoming hadron. The singlet 
part, which involves the gluon distribution function, is formally non-leading in our approach, 
but may nevertheless be important in the HERA kinematic region, where the gluon density 
is high. Estimates of power corrections to singlet structure functions have become available 
very recently [28]; we hope to apply similar techniques to fragmentation functions in a future 
publication. Meanwhile we concentrate on the non-singlet part, where methods similar to those 
used for e + e _ annihilation can be applied. We find that the predicted leading power corrections 
are proportional to l/Q 2 , as in e + e~ annihilation, but their functional forms are different. The 
hypothesis that power corrections are related to a universal low-energy form for the strong 
coupling implies that their magnitudes are given by a single non-perturbative parameter. We 
give quantitative estimates based on the value of this parameter derived from DIS structure 
function data [9]. 

In the following Section we review the approach of Ref. [8]. Sect. 3 presents the standard 
leading-order perturbative treatment of DIS fragmentation, which we modify in Sect. 4 to esti- 
mate the non-singlet power-suppressed corrections using the method outlined in Sect. 2. Some 
numerical results and conclusions are presented in Sect. 5. 



1 



2 Dispersive estimation of power corrections 



We assume that the QCD coupling a s (k 2 ) can be defined down to arbitrarily low values of the 
scale k 2 and that it has reasonable analytic properties, i.e. no singularities other than a cut 
along the negative real axis. It follows that one can write the formal dispersion relation 

*<*•>— jT^*^ (2 ' 1) 

where the 'spectral function' p s represents the discontinuity across the cut, 

Ps{p 2 ) = ^Disc {a s {-p 2 )} = ^ {«s (p 2 e™) - a s (/e~™ ) } • (2-2) 

We now consider the calculation of some observable F in an "improved one-loop" approxi- 
mation, i.e. taking into account one-gluon contributions plus those higher-order terms that lead 
to the running of a s . As discussed in Ref. [8], we expect that 

F = a a (0)?(0) + [°°^4 Ps (p 2 ) HS/Q 2 ) (2-3) 
Jo H 

where the characteristic function F(p 2 /Q 2 ) is obtained by one-loop evaluation of F (divided by 
a s ) with the gluon mass set equal to p [3, 15]. The first term on the right-hand side represents 
the contributions in which a single gluon is produced or exchanged, while the second represents 
those with more complex final or virtual states, e.g. the 'decay products' of a virtual gluon, 
which contribute to the running of a s . In contributions that involve real multi-parton final 
states, Eq. (2.3) with the full spectral function p s in the integrand is obtained only if one sums 
inclusively over a sufficiently wide class of final states. This point will be discussed more fully 
in Sect. 4. 

We can eliminate a s (0) from Eq. (2.3) by means of the dispersion relation (2.1): 

dp 2 
V 2 

Non-perturbative effects at long distances are expected to give rise to a modification in the 
strong coupling at low scales, 5a s , which generates a corresponding modification in the spectral 
function via Eq. (2.2): 

5 Ps (fi 2 ) = ^-Disc {5u s (- p 2 )} . (2.5) 

Inserting this in Eq. (2.4) and rotating the integration contour separately in the two terms 
of the discontinuity, we obtain the following non-perturbative contribution to the observable F: 



r^ Ps (p 2 ) \hp 2 /q 2 ) - m] ■ (2.4) 

Jo P L J 



SF --■ f C °^ r 5a s (p 2 )g(p 2 /Q 2 ) (2.0) 
o p 2 



where, setting p /Q = e 



G{e) = -^-D\sc{H-e)} . (2.7) 

Since 5a s (p 2 ) is limited to low values of p 2 , the asymptotic behaviour of SF at large Q 2 is 
controlled by the behaviour of J-{e) as e — > 0. We see from Eq. (2.7) that no terms analytic at 
e = can contribute to SF. On the other hand for a square-root behaviour at small e, 

T^a^e =► 5F = -^, (2.8) 
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while 

^~a 2 ^elne <5F = a 2 ^| (2.9) 

where 

A = ^ r^M 9 *a.(M 2 )- (2-10) 

Notice that we express the result (2.6) directly in terms of the modification to the strong 
coupling itself, rather than that in the derived quantity a c s which was used in some previous 
publications [8, 9, 20]: 

a ^ = ^^/dV ^ = as{ " 2) - T {^) 2 a ^ 2) + • • • • (2 - n) 

Although a s and a e ff are similar in the perturbative region, they differ substantially at low 
scales, and the former probably has a simpler behaviour. For example, the even moments of the 
effective coupling modification, 

C F dfi_ 2 

have to vanish for all integer values of p, whereas those of 5a s do not. The translation dictionary 
for the moments is in any case rather simple: 

A 2p+ i = (-l) p (p + ±)irA 2p+1 , A 2p = (-l) p pA' 2p , (2.13) 

where 

d . C F f°° dp 2 



A 2p = ^f- % V 2p 6aMM 2 ) , (2-12) 
2-k Jo II 1 



A' 2p = f p A 2p = ^ jf ^ ln^ 2 Sa cS (, 2 ) . (2.14) 

Studies of power corrections to DIS structure functions suggest that A 2 = —A' 2 ~ 0.2 GeV 2 [9]. 

As a clearer representation of the magnitudes of power corrections, we may adopt the ap- 
proach of Refs. [16, 20] and express them directly in terms of moments of a s over the infrared 
region. We substitute for Sa s in Eq. (2.6) 

tfa^-a^-afV), (2-15) 

where a^ T represents the expression for a s corresponding to the part already included in the 
perturbative prediction. As discussed in Ref. [16], if the perturbative calculation is carried out 
to second order in the MS renormalization scheme, with renormalization scale //^, then we have 

afV) = a.(/4) + IbHvl/v 2 ) + k] a 2 ^ 2 R ) (2.16) 

where for Nf active flavours (Ca = 3) 

U€a ~ 2Nf (67-3^)^-10^ 

12vr ' 36vr { ' ' 

The constant k comes from a change of scheme from MS to the more physical scheme [29] in 
which a s is preferably defined at low scales. Then above some infrared matching scale we 
assume that a s (fi 2 ) and aP T (p?) approximately coincide, so that 

m - % r ^ ( as( ^ 2) ■ as{ ^ ] ~ + ^ 

= 2^ ~ as ^ ~ I 61n (^M 2 ) + k + b \ o? a (l4)) , (2.18) 
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where 



(2.19) 



The dependence of <5i on ^ is partially compensated by the /^-dependence of the other terms on 
the right-hand side of Eq. (2.18). The dependence on the renormalization scale //| should help 
to compensate the scale dependence of the perturbative part. Notice that if we take oc Q 2 
then A2 has a logarithmic dependence on Q 2 . In general we do expect power corrections to 
have additional logarithmic Q 2 -dependences (anomalous dimensions), but these are probably 
not given reliably by our 'improved one-loop' approximation. 

3 Fragmentation in DIS 

A complication in DIS, absent from e + e~ annihilation, is the presence in the final state of the 
remnant of the initial-state hadron, i.e. the constituents that did not participate in the hard 
scattering of the lepton. It is expected that the fragmentation of the remnant will be dominated 
by soft, non-perturbative physics. While of interest for studying the hadronization process, the 
remnant fragmentation is not so useful for QCD test, and therefore we concentrate here on 
aspects of fragmentation that are not sensitive to it. This is conveniently done by looking at 
the final state in the Breit frame of reference [30, 31, 32]. 

We consider the deep inelastic scattering of a lepton of momentum I from a nucleon of 
momentum P, with momentum transfer q. The main kinematic variables are Q 2 = -q 2 , the 
Bjorken variable x = Q 2 /2P ■ q and y = P ■ q/P ■ I ~ Q 2 /xs, s being the total cm. energy 
squared. Then the Breit frame is the rest-frame of 2xP + q. In this frame the momentum 
transfer q is purely spacelike, and we choose to align it along the +z axis: 



To a good approximation, the fragmentation products of the remnant will be moving in 
directions close to that of the incoming nucleon, i.e. they will remain in the 'remnant hemisphere' 
H r {p z < 0). On the other hand the products of the hard lepton scattering will tend to be found 
in the 'current hemisphere' H c (p z > 0). In fact in the parton model the scattered parton moves 
along the current (+z) axis with momentum xP + q = ^Q(l, 0, 0, 1). Thus in the parton model 
the current hemisphere looks like one hemisphere of the final state in e + e~ annihilation at centre- 
of-mass energy Q. Fragmentation studies have shown that this similarity is indeed manifest in 
hadron spectra and multiplicities [1]. This makes it natural to define the fragmentation functions 
in terms of particles h appearing in the current hemisphere only, h 6 H c . 

Since we wish to include only particles in the current hemisphere H c , we define the frag- 
mentation function F h for a given hadron species h as a function of the variable z = 2ph ■ q/q 2 , 
which measures the fraction of the hadron's momentum along the current direction and takes 
values < z < 1 in H c : 



P = \Q{l/x, 0, 0, -l/x) , q = ±Q(0, 0, 0, 2) . 



(3.1) 




(3.2) 



The denominator of this expression is the fully inclusive deep inelastic cross section, 



dxdQ 2 



2ttq 2 



{ [l + (1 - y) 2 ] F T (x) + 2(1 - y)F L (x)} 



(3.3) 
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where Ft(x) = 2F\{x) and Fl(x) = Fz{x)/x — 2F\(x) are the transverse and longitudinal 
structure functions (which also have a weak Q 2 -dependence that we do not show explicitly). 
For simplicity we neglect here any contribution from weak interactions (Z° or W ± exchange). 
We shall comment on the effect of this in Sect. 5. The numerator in Eq. (3.2) is the single-hadron 
inclusive cross section, 



d 3 a h 2na 2 



{ [l + (1 - y) 2 ] F%{x, z) + 2(1 - y)F£(x, z)} 



dxdQ 2 dz Q 4 

where Fj. and F^ are generalized structure functions. 
In the parton model (order a°), Fl = F^ = and 

F T (x) = J2e 2 q [q(x)+q(x)]^f(x) 
q 

F*(x,z) = ^2e 2 q [q(x)D q (z)+q(x)D q (z)}, 



(3.4) 



(3.5) 



q(x) and q(x) being the quark and antiquark momentum fraction distributions in the target 
nucleon and D q (z) and Dq(z) their fragmentation functions for hadrons of type h. 




h 



Figure 1: Contribution to jet fragmentation in deep inelastic scattering. 

To first order in a 3 , up to two final-state partons can be emitted in the hard lepton-parton 
subprocess, as illustrated in Fig. 1. The momentum of the struck parton is p = xP/£ (x < £ < 1) 
and we define rj = P ■ r/P ■ q (0 < r/ < 1). The partondevel cross section is 

d 3 a 2vra 2 



dxdQ 2 dr] Q 4 



{ [l + (1 - y) 2 ] F T (x, v ) + 2(1 - y)F L (x, r,)} 



(3.6) 



where (for 77 < 1) 



F ^V) = [ f {CfC^vMx/O +q(x/0\ +T R C i , g (Z,r,)g(x/t)} , (3.7) 
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Cf = 4/3, Tr = 1/2, and g(x) is the gluon distribution. The coefficient functions are [31] 

rj 2 + (1 - rf) 2 



77(1-77) 



CT, fl (e,r?) = [^ + (1-0' 
CL, s (e,r?) = 8^(1-0- (3i 
In the Breit frame P and g are given by Eq. (3.1) and we can write 



where 



p = !Q(l/£,0, 0,-1/0 
r = 5Q(^o,^±, 0,2:3) 

k = ±Q(z ,-z ± ,0,z 3 ) (3.9) 



2 = 2 V - 1 + (1 - rj)/Z 

z 3 = 1 - (1 - 77)/^ 

2 = 1-27? + !?/^ 

?3 = l-?7/C 



2± 



= 2^(1 -77X1 -£)/£. (3.10) 



We can distinguish four subregions of phase space, as illustrated in Fig. 2: 

A: both produced parton momenta k,r in the current hemisphere (2:3,2:3 > 0); 

B: only parton momentum r in the current hemisphere (23 > 0, 23 < 0); 

C: only parton momentum k in the current hemisphere (23 < 0, 23 > 0); 

D: no produced parton momenta in the current hemisphere (23, 23 < 0). 

The 0(a s ) contributions to the generalized structure functions in Eq. (3.4) are of the form 

x{C F Ci, q {Z,rj)q{x/S)[e{z 3 )D q {z/z 3 ) + Q(z 3 )D g (z/z 3 )} 
+ C F Ci, q (t,r,)q(x/0[e(z 3 )D q (z/z 3 ) + Q(z 3 )D g (z/z 3 )} 
+ T R C i , g ^, V )g(x/0[e(z 3 )D q (z/z 3 ) + @(z 3 )D q (z/z 3 )}} (3.11) 

where D g is the gluon fragmentation function. In Eq. (3.11) the coefficient functions Cj )P include 
virtual corrections at 77 = 1, not shown in Eqs. (3.8). 

Changing the variable of integration from 77 to C = 23 or 23 as appropriate in each term, we 
can rewrite this as 



F?(x,z) = fx f'dC 

„ J X J z 



x [K im (£, 0q(x/0D q (z/Q + K i>w (£, Qq{x/£)D g (z/Q 
+ K im (£, Qq{x/0D q {z/Q + K im &, Qq(x/0D g (z/Q 
+ K itgq (Z, Qg(x/0D q (z/Q + K itgq (Z, Qg(x/0D q (z/0] (3.12) 
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where 



Figure 2: Phase space region for jet production in deep inelastic scattering. 



a. 



a. 



a. 



Ki, 9q (t,0 = ^T R c itg (t,z-ZO- 



(3.13) 



It is sometimes useful to convert the convolution integrals in Eq. (3.12) to simple products 
by means of a double Mellin transformation, 

(3.14) 



F(N,M)= f dx ( dzx N - 1 z M - 1 F(x,z) , 
Jo Jo 



so that 



^ k = J2 e l {K im qt> q + Ki, gg qD g + K im qD g + K im qD g + K i>gq gD q + K im gD q ) (3.15) 



where (note the different powers) 



K im (N,M)= f\ f dCt N ( M K im ($,0 
Jo Jo 



(3.16) 



4 Power corrections 



To estimate power corrections according to the method outlined in Sect. 2, we must recalculate 
the relevant observables using a non-zero gluon mass-squared fi 2 = eQ 2 , and then extract 
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the terms that are non-analytic at e = 0. In contributions involving an incoming gluon, the 
"massive" gluon should be treated as an internal line of the singlet process j*q — > q'q'q [28]. 
Thus the cross section involves two massive gluons and is formally beyond the "improved one- 
loop" approximation that we are using. Nevertheless such contributions may be important at 
small x, where the gluon density is large, as is the case at HERA. For the present we avoid these 
complications by considering the non-singlet contribution. For this we require only the quark 
coefficient functions, which become 



CT, q (t,r),e) 



(l-r?)(l -fl + 2fr7(l-7y) 2 -ge 
(1-7,-^e) 2 

(1-^(1-0-^ 



+ 



2^(1 - e) 



4£?7(1 - T]) 2 



in place of those given in Eq. (3.8). The kinematic variables that give the momenta according 
to Eq. (3.9) are now 

z = 2r ? -l + (l-r ? )/£-e 

z 3 = l-(l-r/)/e + e 

z = l-2r] + r]/t + e 

Z3 = l-r]/€-e 

z ± = 2 v /r ? (l-r / )(l-e)/e-er ? . (4.2) 

Thus the phase space region is < 7/ < 1 — e£/(l — £), as illustrated in Fig. 3. The regions 
A,. . . D defined above in terms of the signs of Z3 and Z3 are as indicated. 

The corresponding characteristic functions for the fragmentation functions are given by 
Eq. (3.12) with Ki m (x,z) etc. replaced by (Ci?/27r)/Q i?? (x, z, e) etc., where 

Ki, qq (Z,(,e) = @(l-(-eO@((l-m-C)-eOa, q ((,l-C + a-^,e), 

Kiw&ce) = @(i-e-^)e((i-o(i-e + cc)-< 2 )Q, g (e,e-ec-^,e)- (4.3) 

Note that for brevity we have extracted an overall factor of Cf/2tt, which will be absorbed in 
the non-perturbative factor (2.10). 

Neglecting terms that are analytic or less singular than e In e at e = 0, we find that the Mellin 
transforms (3.16) are given by 

1 1 



}C Tm (N,M,e) 



2S 1 (N) + 2S 1 (M + l)-3 + ± + 1 ^ l 



+ 



+ 



+ 



M + l M + 2 
AS^N + 1) - 45i(M + 1) + 6 + 2N + 2M + 2NM 

M + 2 4 N + A 



In e 



N+ 1 



Kr, qg (N,M,e) 

]C Lm (N,M,e) 
iC L , qg (N, M,e) 



N + 2 
2 



2 

~M + M + l 



M + l 
1 



elne 



M + 2 



lne + 



4 

1+ M 



N + 4 



M + l 



e In e 



N + 2 



e In e 



(4.4) 
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Figure 3: Phase space region with gluon mass-squared /x 2 = eQ 2 . 



where 



N-i 

Si(JV)=E7- 



(4.5) 



Here we have included the virtual contribution, given in Ref. [8]. 

The expression given above for the gluon fragmentation contribution K,T,qg is valid only 
for M > 2. There is an infrared divergence at M = 0, because we integrate the real gluon 
contribution over one hemisphere only, which does not suffice to cancel the divergent virtual 
contribution at ( = 0. For M = 1 there is a contribution of Syfe instead of an elne term, 
implying a l/Q-correction to this moment of the gluon fragmentation function, as is the case 
in e + e~ annihilation [13, 14]. For M = 2 the 1 becomes -1 in the coefficient of elne. All of 
these changes represent extra contributions at the point ( = 0, which we can ignore because the 
fragmentation function at any finite z depends only on the behaviour at £ > z > 0. 

The lne terms in Eqs. (4.4) generate the logarithmic scaling violations in the structure and 
fragmentation functions (see Ref. [8]), while the elne terms give rise to l/Q 2 power corrections, 
as indicated in Eq. (2.9): 




x [H im (Z, Qq{x/i)D q {z/Q + H iiqg (Z, Qq(x/0D g (z/C) 
+ H hqq (tC)q(x/OD q (z/C) + H hqg (tC)q(x/OD g (z/C)] ■ 



(4.6) 
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Inverting the Mellin transforms, Eqs. (4.4) give 



H Tm (£,0 




- 2 - 4£ - +65(1 - + 26' (1 - 5(1 - () 



H L , qq (Z,0 
H T , qg (t,0 



+ [(T^o7~ 4+25 ' (1 ~ c) ] 6(1 ~° 

-<5'(1 - - <5'(1 - + 2<5'(1 - 0<5'(1 - 
[-8^ + 4(5(1-015(1-0 

i-4 + <5(l-0l <5(1 -0+^(1-0- 



(4.7) 



Eqs. (4.6) and (4.7), taken together with the value of A2 — 0.2 GeV 2 and the parton distribu- 
tion and fragmentation functions measured in other processes, provide a quantitative estimate 
of the l/Q 2 corrections to current jet fragmentation in non-singlet DIS. Within the context 
of the dispersive method outlined in Sect. 2, the contributions Hi m from quark fragmentation 
should be reliable, since they are integrated inclusively over all other parton emission. The 
gluon contribution Hx t qg, on the other hand, is less reliably estimated by the "massive gluon" 
approach adopted here. If, for example, a massive virtual gluon splits into two partons, one 
of which goes into the remnant hemisphere and one into the current hemisphere, then only 
the latter will be counted and the full spectral function p s will not be built up in Eq. (2.3). 
The limitations of the massive gluon approach have been studied for l/Q corrections to event 
shapes in full two-loop order [21, 23], and for l/Q 2 corrections to e + e~ fragmentation functions 
in the large- Nf limit, i.e. including only quark loops [14]. In the former case a significant but 
universal enhancement of the power correction was found when going beyond the massive gluon 
approximation, while in the latter case small corrections were obtained, which could however 
become larger when gluon loops are included. In both cases the massive gluon estimate of gluon 
fragmentation effects provided a useful first approximation, and we include it here in the same 
spirit. 

5 Results and conclusions 

To obtain some indicative numerical results from Eqs. (4.6) and (4.7), we assume for simplicity 
that the quark fragmentation function D q is independent of quark flavour and that D q = Dq. 
This is reasonable if heavy flavour production is negligible and one sums over fragmentation 
into all (charged) particles. Then, taking into account the parton-model and non-perturbative 
contributions, we have from Eqs. (3.2)-(3.5) and (4.6) 



where fix) is the charge-weighted parton distribution given in Eq. (3.5) and i?r,g(£) is the 
higher-twist coefficient function for the transverse structure function [8, 9, 10] 1 




-H T , q (0S(l - Q\D q {z/Q + H T , qg {£, ()D g (z/()} 



(5.1) 



HtAO = t-^— - 2 - 4£ + 4(5(1 - + «5'(1 - . 



(5.2) 



Note that the definition here differs from that in Refs. [8, 9, 10] by a factor of — 1/£. 
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It should be noted that the integral of HT, qq (^, C) with respect to £, over the range < ( < 1, is 
not equal to iiy ;? (£), because the region ( < makes a contribution of 25(1 — £) to the latter. 
On the other hand the longitudinal contributions do satisfy the relation 



H Lm (^Q=H L , q (OS(l-0 



(5.3) 



and therefore no longitudinal higher-twist terms appear in Eq. (3.2). Inserting Eqs. (4.7) and 
(5.2) into Eq. (5.1), we find 



F h (z;x,Q 2 ) =D q (z) 



1 v ^i R ^ x - 



(5.4) 



where 



H(z;x) = 2xz 



f(x)D' q (z) J{ X ) ^(z) | D g {z) | 2 



f(x) D q (z) 
4 f 1 

dc 



x 



D q (z) 
/'(*) 



m 
i 



(i-0- 



+x 



+z 



f(x) D q (z) 
D' q (z) 



Jz 

i r 1 



D q (z) 



D q (z) 



X 

? 



'D q (z)f(x)J x — \U ' (5 " 5) 

Figure 4 shows the resulting form of H( Z \ x\ clS ct function of z for various values of x. We 
use the ALEPH [33] parametrizations of the light quark and gluon fragmentation functions for 
charged hadrons at Q = 22 GeV, and the corresponding MRST (central gluon) [34] parton 
distributions. Thus the predictions are at Q 2 = 484 GeV 2 , but H(z;x) depends only weakly 
(logarithmically) on Q 2 , and in any case our method is not reliable at the level of logarithmic 
variations. Results become insensitive to x below the values shown in Fig. 4. Recall, however, 
that we have not computed the singlet contribution, which may well be important at low x 
because of the increase in the gluon distribution there. 

The predicted power corrections are qualitatively similar to those for fragmentation functions 
in e + e~ annihilation [13], though somewhat larger in magnitude. Part of the increase comes from 
the negative higher-twist correction to the transverse structure function in the denominator of 
Eq. (3.2). The rapid rise in the quark contribution at large x and/or z comes from the product of 
derivatives in the first term of Eq. (5.5). The contribution from gluon fragmentation, although 
subject to further corrections as discussed earlier, is estimated to be relatively small for z > 0.2. 

Finally, we comment on the effect of taking weak interactions into account. This introduces 
contributions from the parity- violating structure function Fz(x) in Eq. (3.3) and its analog 
F^(x,z) in Eq. (3.4). However, it was shown in Ref. [9] that the predicted higher-twist contri- 
butions to F3 and Ft are the same, i.e. H^ q = Ht :Q as given in Eq. (5.2). Similarly we find that 



3,99 



Hr,qq and H 3>qg 



HT, q g as given in Eq. (4.7). Hence Eq. (5.5) remains valid, provided 
the parton distribution f(x) is redefined with the appropriate electroweak coefficients in place 
of the charges in Eq. (3.5). 
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Figure 4: Predicted coefficient of A2/Q 2 for fragmentation function in non-singlet DIS. Dashed, dot- 
dashed and solid curves are quark, gluon and total fragmentation. The two sets of curves are for x = 0.3 
(upper) and 0.1 (lower). 
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